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ABSTRACT 

The first extragalactic X-ray binary, LMC X-1, was discovered in 1969. In 
the 1980s, its compact primary was estabhshed as the fourth dynamical black- 
hole candidate. Recently, we published accurate values for the mass of the black 
hole and the orbital inclination angle of the binary system. Building on these 
results, we have analyzed 53 X-ray spectra obtained by RXTE and, using a 
selected sample of 18 of these spectra, we have determined the dimensionless spin 
parameter of the black hole to be a* = 0.92l^QQy. This result takes into account 
all sources of observational and model-parameter uncertainties. The standard 
deviation around the mean value of a=K for these 18 X-ray spectra, which were 
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obtained over a span of several years, is only Aa* = 0.02. When we consider 
our complete sample of 53 RXTE spectra, we find a somewhat higher value of 
the spin parameter and a larger standard deviation. Finally, we show that our 
results based on RXTE data are confirmed by our analyses of selected X-ray 
spectra obtained by the XMM-Newton, BeppoSAX and Ginga missions. 

Subject headings: Galaxies:individual (LMC) - X-rays:binaries - black hole physics 
- binaries:individual (LMC X-1) 



Introduction 



LMC X-1 w as the first extragalactic X-ray binary to be discovered (IMark et al.lll969 
Price et al.lll97ll ). Its X-ray flux is quite constant, varying in intensity by < 25% during 
12 years of moni toring by the All- Sky Monitor (ASM) on board the Rossi X-ray Timing 
Explorer {RXTE; lOrosz et al.l 120091 ). The source ha s been observed by essentially all X -ray 
astronomy missions from Uhuru to Chandra (e.g., iLeong et al.lll97ll : ICui et al.ll2002l ). Its 
bolometric, isotropic luminosity is Lboi ~ 2.2 x 10^^ erg s~^, which is ~ 16% of its Eddington 
luminosity (§5.1). 



A dynamica . 
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Orosz et al. 


2009) 



black-hole model for the system was presented in pioneering work by 



The results of most relevance to this paper are our deter- 
minations of the black hole mass M = 10.91 ± 1.54 Mq and the orbital inclination angle 
i = 36?38 ± 2?02 degree. The distance is likewise important here; we ado pt a distance mod- 
ulus of 18.41 ± 0.10 mag, which corresponds to D = 48.10 ± 2.22 kpc (lOrosz et al.l 120071 . 
20091 ). Accurate values of M, i and D such as these are crucial for the determination of 



black hole spi n by the X-ray continuum-fitting method that we employ (jZhang et al.lll997 
Li et al.ll2005l ). 



To date, our group has published spin estimates for four stellar-mass black holes using 
the X-ray continu um-fitting method : GRO J1655-40, a* = 0.65 - 0.75, and 4U 1543-47 
= 0.75 - 0.85 Jshafee et al.ll2006h : GRS 1915+105 , a* = 0.98 - 1.0 jMcChntock et al 



20061 ): and M33 X-7, = 0.77 ± 0.05 (iLiu et al.ll2008l ). The dimensionless spin parameter 
= cJ/GM"^, where M and J are the mass and angular momentum of the black hole; the 
value a* = corresponds to a Schwarzschild black hole and a* = 1 to an extreme Kerr hole. 
In this paper we show that LMC X-1 harbors a rapidly spinning black hole. Estimates of 
the spins of st ellar-mass black holes are also being obtained by modeling the profile of the 



Fe K line (see Miller et al.ll2002l : |2004J : l2009l . and references therein). 
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In all of our previous work cited above, we have exclusive ly used thermal dominant (TP) 



spectral data in deriving our estimates of black hole spin (IMcClintock fc RemillardI 12006 



Remillard fc McClintockll2006l ). The definition of the TD state is largely based on a spectral 
analysis that makes use of the venerable multi-teraperature and nonrelatiy istic disk model 
known as DISKBB in XSPEC ( Mitsuda et all 1984 : Makishima et al. 1986) plus the stan 



dard power-law component power law ( iTanaka fc Lewin 



1995 



McClintock fc Remillard 



2006: iRemillard fc McClintockl 120061 ): XSPEC is a widely-used X-ray spectral fitting pack- 
age ( ArnaudI 1996 ). The key feature of the TD state is the strong presence of a soft {kT ~ 1 
keV) blackbody-like component of emission that arises in the innermost region of the accre- 
tion disk. The TD state is defined by three criteria: (1) The "disk fraction," the fraction 
of the total 2-20 keV unabsorbed flux in the thermal component, is > 75%; (2) the rms 
power in the power density spectrum integrated from 0.1-10 Hz is < 0.075: ari d (3) QPOs 
are absent or very weak (see Table [2] and text in IRemillard fc McChntoc^ 120061 ). 



Following the selection of the TD sp ectral data, we determi ne the spin parameter a* 
using our relativistic disk model KERRBB2 (IMcClintock et al.ll2006l ). This model incorporates 
a lookup table for the spectral hardening factor / and returns two fit parameters, the spin 
a^, and the mass accretion rate M. In addition to the thermal disk spectrum, there is also 
present a nonthermal tail component of emission, which contributes a few percent of the 1-10 
keV flux. This tail component is widely beli eved to originate from Compto n upscattering of 
the blackbody seed photons in a hot corona (IMcClintock fc Remillardll2006l ). Source spectra 
obtained in the TD state are minimally corrupted by the uncertain effects of Comptonization, 
and they are ideal for the determination of spin via the continuum-fitting method. 



In earlier work (see especially §4.2 in IMcClintock et al.ll2006l ). we have been hampered 
by the lack of an appropriate model of the tail component of emission, which we have 
modeled as a standard power law (powerlaw in XSPEC) or by using one of the existing 
Comptonization models (e.g., comptt). We found both approaches unsatisfactory because 
the power law diverges from realistic Comptonization models at low energies (see §4.2) while 
the existing Comptonization models have more parameters than can be determined. Very 
recently, we have de veloped a sirnple 2- parameter Comptonization model in order to treat 
this tail component (ISteiner et al.ll2009l : §4.21) . The model, which is named SIMPL, has been 
implemented in XSPEC. It is a convolution model that redirects a fraction of the blackbody 
seed photons into a power law. The benefits of using siMPL to model the tail component are 
that it eliminates the unphysical divergence of a standard power law at low energies while 
introducing just two parameters into the fit. Its two parameters are the scattered fraction 
/sc, the fraction of the seed photons that are scattered into a power-law tail component, 
and a photon index P. An illustration of the difference between the models powerlaw and 
SIMPL and a discussion of the benefits of using the latter model are given in §4.2. 
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In this paper, we focus on 53 RXTE spectra obtained using the Proportional Counter 
Array (PCA). These data are central because they are abundant, and especially because 
they allow one to select data via the rms variability criterion (Remillard & McClintock 
2006), which we show is important for the reliable determination of the spin of LMC X-1. 
LMC X-1 was also observed in single observations using XMM-Newton, BeppoSAX, BBXRT, 
and ASCA, and on six occasions by Ginga. In Appendix |Al we analyze the spectral data 
obtained by these five missions and find values of spin that confirm our principal RXTE 
PCA results. 

This paper is organized as follows. In §2 we discuss the selection of our gross sample 
of 53 RXTE spectra and their reduction. In §3 we present the results of a preliminary, 
nonrelativistic data analysis and select a final sample of 18 spectra for analysis. In §4 we 
discuss our relativistic data analysis with emphasis on our first use of the Comptonization 
model SIMPL. We present our results in §5, a general discussion in §6, and our conclusions 
in 57. 



2. Data Selection and Reduction 

As noted in §1, the data of central i mportance in this paper were obtained using the 



large-area PCA detector onboard RXTE (jSwanklll999l ). We consider 53 of the total of 55 



PCA pointed observations that have been performed during the course of the mission; these 
53 observations are identified and characterized in Table 1. We disregard two observations 
that were obtained very early in the mission (1996 February 10 and 1996 March 8) because 
the detector response is poorly known during that period, and the PCA gain was higher 
than at later times. The RXTE PCA data were downloaded from NASA's High Energy 
Astrophysics Science Archive Research Center (HEASARC). 

In our analysis, we only include pulse-height spectra from PCU-2 because fits to the 
power-law spectrum of the Crab show that this is the best calibrated of the five proportional 
counter units (PCUs) that comprise the PCA. Data reduction tools from HEASOFT version 
5.2 were used to screen the event files and spectra. Data were taken in the "Standard 2 
mode," which provides coverage of the PCA bandpass every 16 s. Data from all Xe gas 
layers of PCU-2 were combined to make the spectra. Background spectra were derived using 
the tool pcabackest and the latest "faint source" background model and were then subtracted 
from the total spectra. Redistribution matrix files and ancillary response files were freshly 
generated individually for each PCU layer and combined into a single response file using the 
tool pcarsp. In fitting each of the pulse- height spectra (§§3,4), we used response files that 
were targeted to the time of each observation of LMC X-1. 
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Following the methods spelled out in §3 in iMcClintock et al.l (120061 ) . for the 53 observa- 
tions of the source we individually (1) added the customary systematic error of 1% to all the 
PCA energy channels; (2) corrected the PCA count rates for dead time; and (3) corrected 
the effective area of the PCA using for each observ ation of LMC X-1 th e proximate and 
dead-time corrected observation of the Crab Nebula (IToor fc Sewardlll974j ). The dead-time 
correction factor for LMC X-1 is ~ 1.021. The effective-area correction factor is ~ 1.076, 
which requires us to apply a typical downward correction to the fluxes of ~ 0.930. The net 
effect of the dead-time and effective-area corrections was to reduce the LMC X-1 fluxes that 
we obtained from the analysis of each of our 53 spectra by ^ 5%. 



3. Preliminary Data Analysis 



All of the data analysis and m odel fitting throughout the paper was peformed using 
XSPEC version 12.4 (jArnaudlll996l ). In our two earlier papers that rely primarily on RXTE 
data to estimate the spins of three black holes (jShafee et al.ll2006l : iMcClintock et al.l 120061 ). 
for each source we performed a preliminary analysis using the nonrelativistic disk model 
DISK BB plus the standard power law powerlaw for three reasons: (1) to select the TD 
data (IRemillard fc McClintockll2006l ). (2) to explore the data using this time-tested model, 
and (3) to compare our results with published results. 

For LMC X-1 we have likewise performed a nonrelativistic analysis using this con- 
ventional model, which in XSP EC notation is expressed as phabs(diskbb-I-powerlaw) 
(IRemillard &: McClintockl l2006l ). where phabs is a widely-used model of low-energy pho- 
toelectric absorption. The spectra were fitted over the energy range 2.5-20.0 keV, and 
the XSPEC command "energies" was invoked to accommodate the broader energy range 
required by the convolution model SIMPL. The hydrogen column density was fixed at 
■ /Vh = 4.6 X 10 ^ ^ cm~^, as determ ined by an analysis of Chandra HETG grating data 
(ICui et al.l 120021 : lOrosz et al.ll2009l ). The PCA fits are quite insensitive to the choice of 



the lower cutoff energy because the low-energy absorption is < 10% above 2.5 keV. We 
obtained good fits to this basic model for all 53 spectra without needing to add any of the 
customary additional components (e.g., a Gaussian line or a broad Fe absorption edge). The 
fit results, along with some observational data (e.g., exposure times and count rates), are 
summarized in Table 1. The two pairs of parameters returned by the fit are the two diskbb 
parameters, the temperature at the inner edge of the accretion disk Tin and a normalization 
constant K (columns 7 & 8), and the two powerlaw parameters, the photon index P and 
a normalization constant (columns 9 & 10). 



Columns 2 & 3 in Table 1 identify the observations by calendar date and the midpoint 



- 6 - 



time of the observation. Columns 4 & 5 give the total exposure time and the total source 
count rate. The rms power in column 6, which is a timing parameter, is of special interest, 
and we return to it after commenting on the contents of the remaining columns in the table. 
The four spectral parameters described in the preceding paragraph are listed in columns 
7-10. The disk fraction (DF) as defined in §1, which characterizes the relative strength of 
th e thermal component, is given in column 11. The state of the source (see Table 2 and text 



m iRemillard fc McClintockll2006l ) is given in column 12, and the reduced chi-square and the 
number of degrees of freedom are given in the final column. 

An inspection of column 12 in Table 1 shows that only 3 out of our 53 spectra are 
classified as TD spectra (Nos. 12, 17 & 18)0. We consider our results for these spectra to be 
of the highest reliability, and we therefore refer to them as our "gold" spectra. However, this 
sample is small, and so we consider a second tier of quality, namely those spectra that fail to 
meet the TD criterion that the disk fraction exceed 75%, but which do meet the TD criterion 
that the rms power is < 0.075 (§1). These are our 15 "silver" spectra, and they correspond 
to the following observations listed in Table 1: Obs. Nos. 1-11 and 13-16. Our adopted 
value of the spin parameter (§5.1) is based on the results for all 18 gold and silver spectra; 
in §6 we show that our spin results for both types of spectra are in excellent agreement. 

Meanwhile, the 35 spectra in the lower part of Table 1 (Nos. 19-53) fail to meet the first 
two TD selection criteria. We refer to these as our bronze spectra, and we do not include 
these data in determining our adopted value of a*. However, in §6 we show that these bronze 
spectra yield a value of a^, that is consistent with that obtained from the gold/silver spectra, 
although this value is less well constrained. Throughout the paper we focus our attention on 
our selected sample of 18 gold/silver spectra, all of which satisfy the temporal TD criterion 
that the rms power is < 0.075. 



4. Data Analysis 

The summary in Table 1 of our preliminary analysis results indicates the challenges faced 
in determining the spin of LMC X-1 via the continuum-fitting method. All but three of the 
spectra fail to meet the most fundamental criterion that defines the TD state, namely that 
the unabsorbed thermal component comprise at least 75% of the 2-20 keV flux (column 11). 
A related problem is the steepness of the power-law component (F > 3; column 9), which 



^QPOs were absent in these spe ctra and all 53 RXTE PGA spectra, although they have been reported for 
LMC X-1 in Ginga observations bv lEbisawa et al.l (jl989r ) and in BeppoSAX observations bv lDal Fiume et al. 
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causes the power-law flux to rival or exceed the thermal flux at low energies , a situation 
that we regard as unphysical for a Comptonization model (ISteiner et al.ll2009l ). Finally, 35 
of the 53 spectra fail to meet a second basic criterion of the TD state: rms continuum power 
< 0.075 (§1; Table 1). 

In this section, we describe in turn our workhorse relativistic disk model KERRBB2 and 
the Comptonization model SIMPL. We then summarize the complete set of parameters for 
the composite model that we employ and close by discussing the necessity of flxing the 
power-law index. Meanwhile, in Appendix [Bl we consider and reject the possibility that the 
spectrum of LMC X-1 is contaminated by flux from the nearby X-ray pulsar PSR B0540-69. 



4.1. The Relativistic Accretion Disk Model: KERRBB2 



Our model of a thin accretion disk around a Kerr black hole (ILi et al.ll2005l ) includes 
all relativistic effects such as frame dragging, Doppler boosting, gravitational redshift and 
light bending. It also includes self- irradiation of the disk ("returning radiation") and the 
effects of limb darkening. The effects of spectral hardening are incorporated into the basic 
model KERRBB via a pair of look-up tables for the spectral hardening factor / correspond- 
ing to two representative values of the viscosity parameter: a = 0.01 and 0.1 (see §5.3). 
The entries in this table were cornputed using a second relativistic disk model bhspec 
(jPavis et al.l l2005l : IPavis fc Hubenvl l2006l ). We refer to the model kerrbb plus this ta- 
ble/subroutine as KERRBB2, which is the accretion- disk model that we use hereafter in this 
paper. The model KERRBB2 has just two flt parameters, namely the black hole spin a* and 
the mass accretion rate M, or equivalently, a* and the Eddingto n-scaled bolometric luminos - 



McClintock et al 



ity, Lboi(a*,M)/^Edd, where L^dd = 1-3 x 10^^{M/Mq) erg s"^ (IShaDiro fc Teukolskv 
For a detailed discussion of KERRBB2 and its application, see §4 in 



1983h . 



fl2006f ). 



Throughout the paper, we present detailed results only for the case a = 0.01. However, 
we have computed all of the same results for the case a = 0.1 as well. As we show in §5.3, 
the choice of the higher value of a results in a slightly lower value of a*. Importantly, our 
final adopted value for the spin parameter a* and its error takes into account the inherent 
uncertainty in a (§§5.3,5.6). 



4.2. The Comptonization Model: SIMPL 



In this paper, we make the first application of the Comptonization model siMPL ( ISteiner et al. 



20091 ) in the determination of the spin of a black hole. The model was introduced in §1 where 
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our motivations for develop ing it and its i mport ance in this present work are discussed. The 



model is fully described in ISteiner et al.l (120091 ) where it is applied to data and its perfor 



mance is shown to compare favorably to such physical models of Comptonization as comptt 
and COMPBB. The model SIMPL (SIMple Power Law) functions as a convolution that con- 
verts a fraction fsc of seed photons into a power law with photon index F. There are two 
implementations in XSPEC: SIMPL- 1, which uses a one-sided Green's function and corre- 
sponds to up-scattering only , and siMPL-2, which uses a classical two-sided Green's function 



( ISunyaev fc Titarchuklll980l ) and corresponds to both up- and down-scattering of photons. 
SIMPL- 1 is a "bare-bones" and computationally-efficient implementation. In our application, 
the two models give nearly identical results for the spin parameter. Generically, throughout 
the paper we refer to both versions of the code as siMPL; however, all of the analysis in this 
paper is done using siMPL-1. 

Because siMPL is a convolution model it ties the Comptonized component directly to the 
energy distribution of the input photons, which in the application at hand is specified by our 
disk model KERRBB2. The most important feature of siMPL is that it produces a power-law 
tail at energies higher than those of the thermal seed photons, but the power law does not 
extend to lower energies, which is a property shared by all physical Comptonization models 
(e.g., comptt). Thus, the crucial difference between SiMPL and the standard power law is 
that SIMPL cuts off in a physically natural way, whereas the standard power law diverges 
at low energy. At the same time, siMPL is an empirical model that, while mimicking the 
mechanics of scattering, makes as few assumptions as possible. In particular, the model has 
only two free parameters, as few as the standard power law. 

Figure 1 starkly illustrates the quite different models one obtains employing POW- 
ERLAW versus SIMPL. Shown are fits to one of our 18 selected spectra using the models 
phabs(diskbb+powerlaw) and phabs(simpL(8)DISKBb). Contrasting behaviors between 
POWERLAW and SIMPL are most clearly revealed at low energies in the unabsorbed models, 
which are shown in the rightward pair of panels. Note how the power-law component rises 
without limit, whereas the SIMPL model cuts off in a physically natural way in the manner 
of all Comptonization models. Because powerlaw produces higher fluxes than siMPL at 
low energies, it steals flux from the thermal component and hardens it, thereby reducing 
its normalization constant. Furthermore, the rising power- law component yields inflated 
estimates of A'^h relative to siMPL and other Comptonization models. Thus, although the 
quahty of fit using either of these two-parameter models is comparable, they can yield quite 
different results for both the parameters returned by the thermal component and for A^h- 
These differences between powerlaw and siMPL are most pronounced when the power law 
is relatively steep, F > 3, as in the case of LMC X-1. 
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When the model SIMPL is convolved with a thermal component, in addition to the 
observed thermal component, one also has available the seed spectrum, i.e., the thermal 
spectrum prior to Comptonization. Thus, there are two possible definitions of the disk frac- 
tion: the classical definition of the disk fraction DF that is based on the observed component 
of thermal fiux (Remillard & McClintock 2006), and a new definition based on the fiux of 
the blackbody seed photons, i.e., the model-dependent fiux that emerges from the surface of 
the disk ingoring the effects of the corona that dresses it. We refer to this disk fraction as 
the "naked disk fraction" or NDF. 

The two disk fractions are simply related via the scattered fraction fsc'- DF = (1 — 
fsc)NDF. The choice of which of the two disk fractions to use depends on the purpose at 
hand. The disk fraction DF is currently used to classify the state of the source (§3), whereas 
the naked disk fraction NDF is useful for assessing energetics (e.g., in relating the X-ray 
power to the thermal energy stored in the corona). Throughout this paper, we have chosen 
to use the naked disk fraction and we compute it for the canonical energy interval 2-20 keV. 



4.3. Description of Model Parameters 

In all of the relativistic data analysis, we use KERRBB2 convolved with siMPL in con- 
junction with PHABS, which models the low-energy absorption: i.e., our XSPEC model is 
phabs(simpl®kerrbb2). We use precisely the same energy range, 2.5-20.0 keV, and fixed 
hydrogen column density, A'^h = 4.6 x 10^^ cm~^, that we used in our preliminary data anal- 
ysis (§3). The number of fit parameters is likewise the same, namely four in total: two for 
KERRBB2, the spin parameter a* and the mass accretion rate M, and two for siMPL, the 
power-law index F and the scattered fraction fsc- The normalization constant of KERRBB2 
is fixed at unity, which is appropriate when M, i and D are held fixed at the values quoted 
in §1. For KERRBB2 we also include the effects of limb darkening (lfiag=l) and returning 
radiation (rfi ag=l), and the torq ue at the inner boundary of the accretion disk is set at 



zero {rj = 0; IShafee et al.l l2008a| ]bl). In the case of the Comptonization model siMPL, we 
set switch UpScOnly=l, i.e., we use the siMPL-1 version of the code, which corresponds to 
up-scattering only (§4.2). 



4.4. Fixing the Power-Law Index F 



We further reduced the number of free fit parameters from four to three by fixing the 
spectral index F. We did this because the power-law component is not well constrained, as 
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we now explain. LMC X-1 is classified as a faint source for the PCA, i.e., its total count rate 
is < 40 counts cm~^ s~^ per PCU (Table 1). The PCA background rate is determined using 
a model, in this case the "faint source background model," which is keyed to th e count rate 



of a charged-particle detector aboard the RXTE spacecraft (jjahoda et al.ll2006l ) 



We find that there are quite significant systematic uncertainties in the background count- 
rate model. This is clearly shown by computing the count rates for the 18 gold/silver PHA 
spectra for the energy range 20-40 keV, where our models predict that the source count rate 
is completely negligible relative to the background rate. We first note that these background 
rates indicate that the model background has been overestimated because for 14 of the 18 
spectra the rates are negative. More telling is the fact that these 18 background rates deviate 
from the expected value of zero count rate on average by 4. 60" (std. dev; N = 18). These 
significant deviations show that the background model provides a poor approximation to the 
true background rate. 

We now consider the net source count rates for the 18 gold/silver PHA spectra for the 
energy range 10-20 keV, which is crucial for the determination of the power-law parameters 
(while being completely unaffected by the ~ 1 keV thermal component). Using the 20-40 
keV photon rates as a proxy for the 10-20 keV background rates, we find for half of the 18 
observations (Nos. 1-3, 7-9 and 15-17 in Table 1) that the background rates are fully one- 
fifth of the source count rates. Subtracting this relatively large and uncertain background 
rate from the total count rate significantly corrupts the source spectrum. Thus, for this 
faint source we conclude that we cannot adequately constr ain the power-law index, and w e 



therefore fix it at the reasonable nominal value of F = 2.5 (IRemillard &: McClintockll2006l ). 
The good news is that our final value of the spin parameter a* is very insensitive to this 
choice of F, as we show in 55.2. 



5. Results 

We first present our featured results for the 18 gold/silver spectra in §5.1[ In §5.21 these 
results are compared to those obtained using the standard power law in place of SIMPL and 
also to those obtained using different fixed values of F. In §5.3 we examine the sensitivity of 
our results to varying the metallicity Z and the viscosity parameter a. In §5.41 we consider 
all 53 gold/silver /bronze spectra and examine the correlations between selected pairs of the 
fit parameters. Then, in §5.51 we show how the scattered fraction and the naked disk fraction 
depend systematically on the rms power, and we explore the usefulness of the scattered 
fraction /sc in selecting data. Finally, in §5.61 we derive the error in the spin parameter 
a* that includes all sources of observational uncertainty, as well as the uncertainties in the 
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model parameters. In §§5.1-5.5, the parameters M, i, and D are held fixed at their best 
central values (see §1), and they are allowed to vary only in §5.6. 



5.1. Primary Results 



Our featured result for the spin parameter is given in the bottom line of Table 2: 
= 0.938 ± 0.020 (std. dev.; = 18). Note that this value, which is based on all 18 
gold/silver spectra, is entirely consistent with the value obtained by analyzing just the three 
gold TD spectra (a^, = 0.928 ± 0.014; see footnote c in Table 2). 

Table 2 contains all of our fit results for the 18 gold/silver spectra obtained using our 
relativistic disk model convolved with our Comptonization model. The uncertainties for 
the mean values given in the bottom line are the standard deviation for = 18. From 
left to right, a^, and M are the fitted values of the spin parameter and the mass accretion 
rate, and / is the interpolated value of the spectral hardening factor (§4.1). The two fit 
parameters of SIMPL follow: the photon index F (which is held fixed at 2.5) and the scattered 
fraction /sc (which is fitted). The next quantity is the naked disk fraction as defined for our 
Comptonization model (see §4.2). In the following column, the value of xl the number 
of degrees of freedom show that the fits are all quite good with xl < 1-4. The final column 
gives the Eddington-scaled bolometric luminosity (§4.1) corresponding to the naked disk, 
which is larger than the observed disk luminosity by the factor 1/(1 — /sc), i-e., by 3-5%. 
For all the results given in Table 2 the photon index was fixed at F = 2.5 (§4.4), and the 
column density was fixed at A^h = 4.6 x 10^^ cm~^ (§3). 

The luminosity is a modest fraction of the Eddington value and shows little variation over 
the 1996-2004 period during which these 18 spectra were obtained: L/L-^dd = 0.160 ± 0.007 
(Table 2). The relative constancy of LMC X-1 is also attested to by the 12- year record of 
the RXTE All-Sky Monitor: The intensity of LMC X-1 has been more or less steady at 
K, 20 mCrab (1.5- 12 keV) with variations averaged over 10-day intervals of only f« ±25% 
(jOrosz et al.ll2009l ). In this work, as in all previous work, we select only those observations 
for which L/L-^dd < 0.3, which corresponds to a disk that is geomet rically thin at all radii, 
with a height-to-radius ratio less than 0.1 (IMcClintock et al.ll2006l ). An inspection of the 
last column in Table 2 shows that all 18 gold/silver spectra easily meet this requirement. 

On a time scale of many thousands of years, one fully expects the spin of the black 
hole to be constant, essentially u naffected by either accretion torqu es or any process that 
might extract angular momentum (IKing &: Kolblll999l : lLiu et al.ll2008l ). As shown in Figure [2] 
and as expected, the spin of LMC X-1 over the 8-year span of the RXTE observations is 
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consistent with a constant value of a* to within ^ 2% (std. dev.; N=18). 

Importantly, the high spin we find using RXTE PCA data, a* = 0.938 ± 0.020, is 
confirmed by X-ray missions with low-energy response that reaches well below the 2.5 keV 
cutoff of the PCA. As we show in Appendix A, when we apply the data selection criterion 
developed in §5.5 (/sc < 4%), we find values of spin using spectra obtained by the XMM- 
Newton, BeppoSAX, and Ginga missions that are consistent with those we obtain using our 
RXTE spectra. 



5.2. Results Obtained Using Other Models and Using Different Values of F 

In all earher work, we have relied on the standard model phabs(kerrbb2-(-powerlaw) 
in deriving th e spin parameter (^3). We found that this model was generally satisfactory 



(but see §4 in iMcClintock et al.l l2006l ) because the power law was not too steep (F ~ 2.5). 
We now apply this old model to the 18 gold/silver spectra of LMC X-1 with F > 3 (Table 
1) in order to see how much the steepness of the power law affects our determination of 
a*. In making this comparison, we follow exactly the steps in our earlier analysis except 
that we model the Comptonized tail component of emission using the additive component 
POWERLAW rather than convolving KERRBB2 with SIMPL. For example, we fix F = 2.5 
and A^H = 4.6 x 10^^ cm~^ and fit over the energy range 2.5-20.0 keV. The resultant mean 
fit parameters for the 18 spectra are given in line 3 of Table 3. For ease of comparison 
our featured mean fit parameters, which are repeated from Table 2, are given in the first 
line of Table 3. Comparing lines 1 & 3, one sees that the old model yields a value of 
a^, = 0.953 ± 0.018, which differs from our featured value of a* = 0.938 ± 0.020 by only 0.7 
standard deviation. 

How sensitive are our adopted results in Table 2 to our chosen value of F = 2.5 (§4.4)? 
In order to answer this question, we have recomputed the results given in Table 2 with the 
photon index fixed at F = 2.3 and F = 2.7, and we find respectively, a* = 0.944 ± 0.019 and 
= 0.932 ±0.021. These values differ from our adopted value of a* = 0.938 ±0.020 by only 
0.2 standard deviations. Even if we fix the photon index at F = 3.1, the mean value in Table 
1 for the 18 gold/silver spectra (as well as for all 53 spectra), we find a* = 0.918 ± 0.024, 
which differs from our adopted value by 0.6 standard deviation. This small difference is 
not unexpected given that so few thermal photons are scattered into the tail component, 
/sc = 0.036±0.005 (Table 2). That a* is anti-correlated with F is expected because a steeper 
power-law component steals photons from the high-energy tail of the thermal component, 
and the resultant softening of this component implies a lower value of spin (and vice versa). 
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As a final detail, we note that for just 4 of our 18 spectra there is clear evidence for an 
Fe emission line in the fit residuals. In deriving our primary results presented in §5.1, we did 
not include an Fe line because the fits are quite good without it. As a test, we have refitted 
all 18 spectra while including the Fe line, and we find a resulting mean spin parameter of 
a* = 0.922 ± 0.021, which is quite consistent with our featured result above (see footnote e 
in Table 2). 



5.3. Results for Different Values of Metallicity and Viscosity Parameter 



^The metallicity of the LMC, and presumably LMC X-1 as well, is Z ^ 0.3^0 (jPopita fc Sutherland 

20031 ) . whereas our results are all given for Z = Zq. How important is this difference in 



metallicity? We recalculated all of the results summarized in Table 2 for the 18 gold/silver 
spectra using readily-available BHSPEC table models that were computed for Z = O.IZq. 
We found = 0.937±0.020 (std. dev.; N = 18), which is almost identical to the bottom-line 
result quoted in Table 2 {a^ = 0.938 ±0.020). Thus, varying the metallicity has a completely 
negligible effect on our results. 

In order to streamline the paper, throughout we quote detailed results only for one 
value of the viscosity parameter, a = 0.01. Presen tly, there is an inherent uncertainty 
in a, and we direct the reader to iKing et al.l (120071 ) for a more comprehensive review of 



numerical and observational constraints. I n addition, a detaile d discussion of BHSPEC's 
spectral dependence on a can be found in iDone fc David (120081 ). A low value is motivated 
by shearing box sim ulations with zero net magnetic flux, which show a approaching zero as 
resolution increases (jPessah et al.l 120071 : iFromang fc Papaloizoul 120071 ). The spectral models 
are remarkably insensitive to a for a < 0.01. Therefore, results obtained for models with 
a = 0.01 are representative of those obtained for models with even lower values of a. 

Meanwhile, larger values of the viscosity parameter, 0.1 , are suggested by global 



GRMH D simulation s (IHawlev fc KrolikI 12001 



Shafee et al.l l2008a| ) and the disk instabihty 



model (lLasotall200ll ). Our results for a* are modestly sensitive to these larger values of a. 
Therefore, we have recalculated all of the results summarized in Table 2 using BHSPEC 
table models computed for a = 0.1. In this case, we find a slightly lower value of the spin 
parameter: a^, = 0.908 ± 0.023. This value is 1.0 standard deviations less than the value 
given in the bottom line of Table 2 for a = 0.01 (a* = 0.938 ± 0.020). In deriving our 
final adopted value of the spin parameter and its error (§5.6), we take into account the 
uncertainty in the viscosity parameter by computing our results using both of our fiducial 
values of a: 0.01 and 0.1. 
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5.4. Graphical Presentation of Results for All 53 Spectra 

Our adopted results (§5.1) are based exclusively on the 18 gold/silver spectra. All of 
the fitting results for these spectra are given in Table 2, and histograms for the four fit 
parameters are shown as filled blocks in Figure [31 Although we do not present a table of our 
results for the 35 bronze spectra (Nos. 19-53 in Table 1), we do show histograms based on 
the results of our analysis of these spectra as open blocks in Figure [3l which are superposed 
on the solid-block histograms. An inspection of the figure shows that, when compared to 
the 18 gold/silver spectra, the fits to all 53 spectra yield on average higher values of the spin 
parameter a^,, lower mass accretion rates (M), more scattering of the thermal photons (/sc) 
and smaller naked disk fractions (NDFs). A quantitative comparison of these differences is 
provided by a comparison of the results for the 18 gold/silver spectra, which are given in the 
top line of Table 3, with the comparable results we obtained for all 53 spectra, which are 
given in the bottom line of Table 3. 

Figure m shows correlation plots between pairs of fit parameters. The 18 gold/silver 
spectra are indicated by the filled circles and the 35 bronze spectra by open circles. Figure 
shows a strong correlation between a* and the mass accretion rate M. The range of the 
correlation is much larger for the bronze spectra, extending to a* ~ 1. This correlation is 
rather uninteresting because it simply reflects the fact that M is a normalization constant 
that must necessarily decrease if the measured parameter a^, fluctuates upward (or vice 
versa). In Figure© we recast this correlation by plotting a* versus rjM, where the efficiency 
parameter rj = Lboi/Afc^ is computed for the particular value of a*. As shown, there is no 
correlation between and r/M, i.e., between a^, and the bolometric luminosity. 

Figure shows a tight correlation between the naked disk fraction and the scattered 
fraction. Our adopted disk fraction (§4.2) is the ratio of the flux in the naked disk component 
-Fnd to the total unabsorbed flux from 2-20 keV, which is the observed thermal flux plus 
the power-law component, (1 — /sc)-^nd+-Fpl- As /sc increases and thermal photons are 
directed into the power-law tail, FpL is expected to increase substantially because each 
Compton-upscattered photon is boosted in energy. Thus, as /sc and FpL increase, the disk 
fraction must decrease, as shown. 

The relationship between a* and the scattered fraction is shown in Figure Ad. There is 
no clear correlation for the 18 gold/silver spectra. However, among the 35 bronze spectra 
is seen to increase with scattered fraction, although the correlation is weak and there are 
two discrepantly low points. 
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5.5. Using the Scattered Fraction to Select Data 

Figure [5] shows correlation plots of the scattered fraction and naked disk fraction versus 
the rms power (§1). The solid vertical line at rms power = 0.075 divides the bronze data 
on the right from the gold/silver data on the left. These latter data are weakly correlated, 
showing a tendency for the scattered fraction to increase and the naked disk fraction to 
decrease as the rms power increases. To the right of the dividing line the sense of the 
correlations are the same, but the correlations steepen and become very pronounced. Of 
course, one expects the scattered fraction and the naked disk fraction to correlate in the 
opposite sense given the tight anti-correlation between these quantities that is shown in 
Figure Hb. 

The TD-state criterion that the rms power is < 0.075 has played a central role in the 
selection of our data (§3). Unfortunately, the use of the rms power criterion, which requires 
timing data up to ~ 10 Hz and is based on the particular broad energy response of the 
PCA, is applicable to RXTE data only. Meanwhile, the good correlation between the rms 
power and the scattered fraction (Figure [5]) allows the latter quantity to be used in selecting 
data as a proxy for the rms power. (One could use the naked disk fraction instead of /gc; 
however, we favor /sc because it is a readily-available fit parameter.) We now investigate 
this possibility in the case of LMC X-1 using the model PHABs(siMPL(g)KERRBB2). 

In the top panel of Figure O we show the rms power for all 53 spectra plotted in order of 
increasing rms power. The horizontal line marks the TD-state criterion rms < 0.075 and the 
vertical line separates our 18 gold/silver spectra on the left from the 35 bronze spectra on the 
right. With the same ordering of the spectra, we show in the lower panel for all 53 spectra 
the scattered fraction. The dashed line for /gc = 4% was determined as follows: We assumed 
that both the rms power and /gc follow normal distributions and are uncorrelated with the 
other fitted parameters. For the rms power we find a mean value of rms = 0.081 ± 0.016 
(std. dev.; n = 53), which is 0.38(7 above the rms = 0.075 line in the top panel. Given the 
correlation between fsc and rms power (Figure E]), we expect the corresponding scattered- 
fraction cut to also be 0.38a below the mean value, i.e., at f^c = 4% (rounded to the nearest 
percent); this is the level indicated by the dashed line in the lower panel. 

The solid symbols in the top panel correspond to the 16 spectra in the lower panel that 
have fsc < 4%. It can be seen that all of the spectra selected by scattered fraction are 
consistent or nearly consistent with the TD criterion rms < 0.075. Likewise, in the lower 
panel one sees that most (14 out of 18) of the gold/silver spectra have fsc < 4% while 4 of 
them have 4% < fsc ^ 5%. Meanwhile, most of the 35 bronze spectra with rms significantly 
greater than 0.075 are soundly rejected by the criterion fsc < 4%. 
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We now compare the fitting results obtained for tlie data selected using the criterion 
rms < 0.075 to that selected using the criterion fsc < 4% and find that they are essentially 
identical: The spin parameter and the scattered fraction for the rms-selected data are re- 
spectively a* = 0.938 ±0.020 and U = 0.036 ±0.005 (Table 2), and for the /sc-selected data 
are a, = 0.938 ± 0.016 and fsc = 0.033 ± 0.003. These results are based only on LMC X-1. 
However, they indicate that whenever high-quality timing data is unavailable for a source 
that it may be possible to establish a useful scattered-fraction criterion capable of rejecting 
data with high flicker noise. Of course, the scattered-fraction cut level that one adopts may 
vary from source to source. 

In Appendix A, we make the first use of our newly-defined scattered fraction criterion, 
fsc < 4%, in order to select data obtained by five X-ray missions for which the RXTE rms 
data-selection criterion cannot be apphed. 



5.6. Comprehensive Error Analysis: Final Determination of Spin Parameter 



The statistical error in the spin parameter is quite small: a^, = 0.938 ± 0.020 (a = 0.01; 
Table 2). At this level of precision, the observational error is dominated by the uncertainties 
in the input p arameters: the bl ack hole mass M, the orbital inclination angle i and the 



distance D (§l: IOrosz et al.ll2009l ). In order to determine the error in a^, due to the combined 



uncertainties in M, i and we performed Monte Carlo simulations assuming that the 
uncertainties in these parameters are normally and independently distributed. Specifically, 
we (1) generated 3000 parameter sets for M, i and D] (2) computed for each set the look- 
up table for the spectral hardening factor / using the model bhspec; and (3) using these 
/-tables, we obtained by fitting our standard model (§4.3) to the 18 gold/silver spectra. 

The results for the 3000 simulations are shown in Figure [3 Panels a-c show the dis- 
tribution of the values of a^, for the parameters M, i and respectively. Each set of 18 
spectra is averaged and represe nted by a single d ata point. In each panel, the value of the 
parameter in question is given (jOrosz et al.ll2009l ). and a large central dot is plotted among 
the 3000 distributed points at a* = 0.937, which is the median value of the distribution. The 
histogram for the displacements in a=K about the most probable value is shown in panel d. As 
indicated, the 1 a error in a^, is (±0.040, -0.061), which is much larger than the statistical 
error of ±0.020 (Table 2). The error is dominated by the uncertainty in M; the uncertainties 
in i and D are relatively unimportant. 



The analysis above is for our standard value of a = 0.01. We also performed the Monte 
Carlo analysis for our other fiducial value of the viscosity parameter, a = 0.1, which as we saw 
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in §5.3 results in a slightly lower mean value of the spin parameter, a* = 0.908 ±0.023. (The 
computed median value of the distribution is = 0.907.) In Figure [H we show the results 
of this analysis as a dashed-line histogram of spin displacements; the histogram for a = 0.01, 
which is copied from Figure Wd, is also shown by the smaller solid-line histogram. The 
summation of these two histograms results in the large histogram (heavy solid line). This 
combined distribution, which corresponds to a total of 6000 simulations, has a median spin 
value of 0.92 and implies a la error of (+0.05, -0.07). Thus, considering both the uncertainties 
in M, i and D and the inherent uncertainty in the value of the viscosity parameter we arrive 
at our final result for the spin of LMC X-1: a* = 0.92]'^o;q7 (Icr). 



6. Discussion 

Throughout, we have focused on the three gold spectra, which fully satisfy the TD-state 
criteria, and 15 additional silver spectra. These latter spectra meet the TD-state rms power 
criterion. However, they yield low values of the disk fraction (Table 1) that fail to meet 
the TD-state disk-fraction criterion because the standard power-law model is inadequate for 
modeling their steep Compton-tail components. We have shown that these 18 gold/silver 
spectra when analyzed using the Comptonization model SIMPL yield large values of the naked 
disk fraction (Table 2). Restricting ourselves to this sample of 18 spectra and considering only 
the statistical error, we obtained our featured value of the spin parameter: a* = 0.938 ±0.020 
(std. dev.; N=18). For comparison, if we consider all 53 spectra (i.e., including the bronze 
spectra) we obtain: a* = 0.961 ± 0.033 (std. dev.; N=53; Table 3). That is, we obtain 
a statistically consistent result, although with a somewhat higher spin and with a larger 
uncertainty. We have chosen to base our result on just the 18 gold/silver spectra because 
they meet the TD-state rms criterion and because their spectra are strongly dominated by 
the thermal component, i.e., their 2-20 keV naked disk fractions are > 86% when computed 
using SIMPL (Table 2). When we include all the sources of observational error and model 
parameter uncertainties, our final result for the spin parameter is a* = 0.92]'^o;g7. 

The greatest uncertainty in our spin estimate arises from the uncertainties in the va- 
lidity of the disk model we employ. For example, the model assumes a razor-thin disk 
for which the viscous torque is assumed to vanish at the innermost stable circular orbit 
(ISCO ), an assumption that, it has been argued, does not hold in the presence of magnetic 



fields (IKrolik fc Hawleyll2002l ). However, our preliminary hydrodynamic work indicates that 
the effects of these viscous torques and emission from within the ISCO are modest for the 
thin disks (H/R < 0.1) and low luminositi es (L/LEdd < 0.3) that we restrict ourselves to 



flMcClintock et al.l l2006l : IShafee et al.ll2006l ). In the case of LMC X-1, with its modest lu- 
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minosity (L/LEdd ~ 0.16; Table 2) and hig h spin, these hydrod ynamic results imply a very 
small uncertainty, Aa* < 0.01 (see Fig. 9 in lShafee et al.ll2008bl ). Another uncertainty lies in 
our assumption that the black hole spin is approximately aligned with the angular momen- 
tum vector of the binary. As Figure 36 indicates, if any misalignment is < 2°, then it will 
contribute an error in the spin parameter a* that is no larger than our total observational 
error. A future X-ray polarimetry mission may allow a check on our assumption of alignment 
Li et al.ll2009[ ). 



Our preliminary hydrodynamic results are supported by our more recent three-dimensional 
general relativistic M HD simulations of an appropriately thin disk with H/R ~ 0.05 — 0.1 
( IShafee et al.l l2008al ). This model is for a non-spinning black hole. In steady state, the 
specific angular momentum profile of the inflowing magnetized gas deviates by < 2% from 
that of the standard thin disk model of iNovikov fc Thornd ( 119731 ). In addition, the magnetic 
torque at the ISCO is only ~ 2% of the inward flux of angular momentum at this radius. 
These results indicate that magnetic coupling across the ISCO is relatively unimportant for 
geometrically thin disks, at least in the case of a non-spinning black hole. 

In short, our recent theoretical studies indicate that the disk emission is truncated 
rather sharply at the ISCO and that measuring the truncation radius provides a reliable 
estimate of spin. In addition, there is a long history of observational evidence suggesting 
that fitting the X-ray continuum is a promising approach to measuring black hole spin. 
This history begins in the mid-1980s with t he simple non-relativistic multicolor disk model 
( iMitsuda et al.lll984j : iMakishima et al.lll986l ) now known as diskbb, which returns the color 
temperature at the i nner-disk radius Rin (§3). In their review paper on black hole binaries, 
Tanaka fc LewinI (119951 ) summarize examples of the steady decay (by factors of 10-100) of 
the thermal flux of transient sources during which R^^ remains quite constant (see their Fig. 
3.14). They remark th at the consta r icy of i?in suggests that this fit parameter is related to 
the radius of the ISCO. IZhang et al.l (119971 ) then outlined how, using a relativistic disk model 
and corrections for the effects of radiative transfer, the fixed inner disk radius provides an 
observational basis to infer black hole spin. More recently, this evidence for a constant inner 
radius in the thermal state has been presented for a number of sources via plots showing that 
the bolometric luminosity of the thermal compo r ient is approximately proportional to T-^ 



(iKubota et al 



2001 



McClintock et al. 



Kubota fc Makishimal 12004 : ICierhhski fc Pond 12004 : lAbe et al.l 12005 



20091 



We now consider the possibility that our results are affe cted by a component of ab- 
sorption that is not included in our model. ICooke et al.l (120081 ) have found evidence for an 
ionization cone within a few parsecs of LMC X-1. They argue that the cone must be pow- 
ered primarily by photoionzation from the accretion disk, but require a 13-300 eV photon 



- 19 - 



flux which exceeds the extrapolated hard X-ray model by a factor ~ 10. Both this excess 
and the anisotropy of the ionized gas suggest that the X-ray emission might be partially 
obscured near the source, possibly by the outer disk. If this hypothesis is correct, we could 
be underestimating the true disk luminosity and, therefore, overestimating the spin. How- 
ever, we stand by our estimate that the bolometric luminosity of the disk at the source is 
L/L-Edd ~ 0.16 (Table 2) while regarding the evidence that the disk is substantially more 
luminous as tentative for several reasons: e.g., the possibility of accretion rate variations 
over the recombination time of the ionized gas; the lack of evidence for obscuration in the 
X-ray spectrum; the close proximity of other external ionizing sources; and the likelihood 
of significant EUV flux via X-ray reprocessing in the accretion disk. On this last point, 
our quantitative estimate indicates that reprocessing can provide the ionizing photon flux 
required to explain the He I I A4686 line flux: Using the accretion-disk reprocessing model 
described in Appendix A of lOrosz et al.l (120091 ) . and assuming that the annular size of the 
reprocessor is i? = 100 - W^GM /c^, we flnd a 54.4-3 00 eV photon flux of 4 x 10^^ photons 



s ^, which matches the flux that lCooke et al.l (120081 ) flnd is required to power the ionization 



con 



What is the origin of the high spin of LMC X-l's black hole primary? Was the black hole 
born with its present spin, or was it torqued up gradually via the accretion flow supplied 
by its companion? In order to achieve a spin of a* ~ 0.9 via disk accretion, an initia lly 
non-spinning black hole must accrete more than 4Mq from i ts donor (IKing fc Kolblll999l ) in 
becoming the M = 1O.9M0 black hole that we observe today (jOrosz et al. I l2009h . However, to 
transfer this much mass even in the case of Eddington-hmited accretion {MecLcl = LecIci/c^ ~ 
2.50 X 10~^Mo/yr) re quires > 15(3 milli on years, whereas the age of the system is only 
about 5 million years (jOrosz et al.l 120091 ). Thus, we conclude that the spin of LMC X-1 
is nat al, which is the same conclusion that has been reached fo r three other stellar black 
holes jMcChntock et allbood : Ishafee et allbood : Ihiu et al.l[2008h . We note, however, that 
an alternative scenario for the origin of high s pins for black holes in X-ray binaries has been 
proposed that invokes hypercritical accretion (IMoreno Mendez et al.l 120081 ). 



^We obtained the value 4 x 10^^ by integrating the photon flux, assuming blackbody emission, over the 
energy range from 54.4 eV to 300 eV and over the entire radial extent of the reprocess ing region. We use d 
the radial temperature profile specified by Equations (A3) and (A5) in Appendix A of lOrosz et al.l (|2009r ). 
Our result is insensitive to the values we assumed for the inner and outer radii of the annular reprocessor. 
For example, if we halve the inner radius or double the outer radius, the 54.4 eV to 300 eV photon flux 
changes by only « 10%. 
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7. Conclusion 

We find tliat tlie spin of tlie black hole primary in LMC X-1 is a* = 0.92^Q Qy, where the 
uncertainty includes all sources of observational error and model parameter uncertainties. 
This result is based on an analysis of three thermal dominant (TD) spectra, which we refer to 
as our gold spectra, plus 15 silver spectra. The observations that yielded the 15 silver spectra 
satisfy the TD criterion that their rms power is < 0.075. However, when analyzed using the 
standard power-law model they fail to meet the TD-criterion that their disk fractions are 
> 75%. We show that the failure of the silver spectra to meet this criterion is caused by 
the unphysical behavior of the standard power-law model, which diverges at low energies 
and therefore intrudes on the thermal component of emission. We sidestep this problem 
by using our Comptonization model SIMPL. If we admit an additional 35 bronze spectra 
(with rms power > 0.075) into the analysis and consider all 53 RXTE spectra, we obtain a 
consistent, somewhat higher, and less certain value of the spin parameter. Meanwhile, our 
spin results for the 18 gold/silver RXTE spectra are confirmed by our analyses of XMM- 
Newton, Beppo SAX and Ginga spectra (see Appendix A). 
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A. Determinations of the Spin Parameter for Five Additional X-ray Missions 

As stated in the final paragraph of §5.1, our featured result for the 18 gold/silver RXTE 
spectra for the case a = 0.01, a* = 0.938±0.020 (std. dev.; N=18), is confirmed by an analysis 
of data from three other X-ray missions. In the paragraphs below, which are prefaced by 
a mission name, we consider in turn X-ray spectral data for LMC X-1 that were obtained 
during the course of the five missions named in Table H] (column 2). With the exception 
of the BBXRT and Ginga data, all the data were downloaded from the NASA HEASARC. 
As discussed below, after applying the data selection criterion developed in §5.5, namely 
scattered fraction fsc < 4%, we find that only three of these missions provide reliable 
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estimates of spin (Table HI lines 1-3), i.e., estimates of spin that are comparable in quality 
to the 18 gold/silver RXTE PCA spectra. As in the analysis of the RXTE data, for all 
five of these missions we fix the power-law index at 2.5, which is a necessity for four of the 
five missions whose detectors are unresponsive above ~ 10 keV. As indicated below, small 
corrections w ere applied to the effe ctive areas of the various detectors to align them with 
the standard iToor fc SewardI (119741 ) spectrum of the Crab. We consider here only the case 
a = 0.01. 

XMM-Newton: A single 4940 s observation was performed on UT 2000 October 21 
using the EPIC-pn camera in timing mode. The data, which are free of pile-up effects, were 
reduced using precisely the proced ures recommended for timing dataB We fitted the data 
over the energy range 1.0-10.0 keV (j Wilms et al.ll2003l ) and corrected the detector's effective 
area by the factor 0.889 (i.e., the effective a rea is increased by dividing by the factor 0.889), 
which is based on observations of the Crab (IKirsch et al.ll2005l ). Conclusion: The data meet 
the requirement /sc < 4%, and the spin is consistent with the RXTE value to within 0.9 
standard deviations (Table HI line 1). 

BeppoSAX: A single observation was made on UT 1997 October 5 using the narrow- 
field instruments aboard BeppoSAX. The four de tector systems , as w ell as this specific ob- 
servation of LMC X-1, are described in detail in iHaardt et al.l (120011 ). We have chosen to 
analyze only the data obtained by the Medium Energy Concentrator Spectrometer (MECS 
= MECS2 + MECS3) because the PDS and HPGSPC spectra are noisy and the LECS 
spectrum lacks a proper response file. The MECS data have been reduced following the 
standard procedures described in the Cookbook for the BeppoSAX NFI Spectral Analysifl 
The net eposure time was 40 ks. The spectra obtained by the two MECS detectors, MECS2 
and MECS3, were fitted jointly over the energy range 1.65-10.5 keV; no correction to the 
effective area was required. Conclusion: The data meet the requirement /sc < 4%, and the 
spin is consistent with the RXTE value to within 1.5 standard deviations (Table HI line 2). 

Ginga: The Ginga Large Area Proportional Counter (LAC; 1.2-37 keV) observed LMC 
X-1 on six UT dates: 1987 April 22, 1987 July 16, 1987 September 30, 1988 August 15, 
1990 June 21 and 1991 February 15. The respective exposure times were approximately 
7290, 12831, 4841, 315, 326, and 266 s. We extracted the spectra following the procedures 
described in lEbisawa et al.l (119891 ). We applied to each spectrum an effect ive-area correction 
factor of 0.968 based on an observation of the Crab ( Turner et al.lE989l ). We fitted all six 
spectra; however, in Table |4] we report our results only for the spectrum obtained on 1987 



'^ftp://ftp.xray.mpe.mpg.de/xmm/service/data_analysis/EPIC J^N/timing/timing_mode.html 
"'http: / /heasarc. gsfc.nasa.gov/docs/sax/abc/saxabc/saxabc. html 
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September 30, the only spectrum that satisfies our data selection criterion (/sc < 4%). This 
spectrum was fitted over the restricted energy range 1.2-12.0 keV because the spectrum 
contains large and unexplained residual features at higher energies. Conclusion: The data 
meet the requirement /sc < 4%, and the spin is consistent with the RXTE value to within 
0.8 standard deviations (Table HJ line 3). 

BBXRT/Astro-1: Four observations were performed on UT 1990 December 6 & 7. 
We disregard the observations on December 7 because of their high background levels and 
large off-axis angles; we consider only the pair of observations performed on December 6 
with exposu re times of 1235 s an d 616 s. The data have been reduced following the standard 
procedures (jSchlegel et al.lll994l ). The two spectra were fitted jointly o yer the energy rang e 
0.7-7.0 keV, and the effective area was corrected by the factor 0.973 (IWeaver et al.lll995l ). 
Conclusion: The scattered fraction fails to meet the requirement /sc < 4%, the spin value 
is unreliable, and we disregard it (Table HJ line 4). However, we note that the BBXRT spin 
value is consistent with the RXTE Ydlne to within 1.0 standard deviation. 

ASCA: A single 3465 s observation was performed on UT 1995 April 2. We consider 
only the data obtained by the Gas Imaging Spectrometer (CIS). We extracted the spectra 
following the standard procedures described in the ASCA ABC guide 0. The two CIS spectra 
(GIS2 and GIS3) were fitted jointly over the energy range 0.7-8.0 keV. No correction to 
the effective area is re quired because the GIS effective area calibrations were based on the 
Toor &: SewardI (Il974l ) spectrum of the Crab. Conclusion: The scattered fraction fails to 
meet the requirement /sc < 4%, and the spin value is unreliable. We therefore disregard it 
while noting that it is exceptionally low, a result for which we have no explanation. 



B. Contamination by the Pulsar PSR B0540-69 



The pulsar PSR B0540-69 is located 25' from LMC X-1. iHaardt et all (1200 If ) concluded 
that their BeppoSAX PDS spectra (15-100 keV) of LMC X-1 were contaminated by the 
presence of this nearby pulsar. Even though the PCA has a narrower field-of-view than 
the BeppoSAX FDS, 1?0 versus 1?5 (FWHM), our spectra of LMC X-1 could likewise be 
significantly contaminated. We now describe three lines of evidence showing that this is 
not the case: We first estimate the contribution of the pulsar to the total fiux using its 
well-determined spectrum; we then search for the expected signature of the pulsar's hard 
spectrum in our data; and lastly we subtract a faked version of the pulsar's spectrum from the 
source spectrum. (We disregard PSR J0537-6910, which also lies in the PCA field of view. 



■"^http: / /heasarc. gsfc.nasa.gov/docs/asca/abc/abc. html 
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becau se it is further off-axis and significantly fainter than PSR B0540-69; ICusumano et al. 
1998h . 



PSR B0540-69 and its 2" —2)" pleri on nebula were discovered 25 years ago by the Einstein 
X-ray observatory (jSeward et al.lll984l ). The source has since been closely observed by many 
missions (e.g., ROSAT, RXTE, Chandra, and INTEGRAL). It is quite similar to the Crab. 
Like the Crab, all studies to date show that its hard sp ectrum can be fitted with a single 



powe r law with a phot on index F in the range 1.8 to 2.2 (ISunyaev et al.lll99Cll : iKaaret et al. 



200ll : iGotz et al.ll2006l ). On the other hand, t he power-law component in the spectru m of 
LMC X-1, is much softer, T ~ 2.5 - 3 (Table 1; iNowak et allboOll : Iflaardt et ahlboOlh . We 
now estimate at what energy the flux in the fiat pulsar spectrum first exceeds the flux of 
LMC X-1. With LMC X-1 centered in the PCA, the collimator transmission factor for the 
pulsar is T = 0.57. Adopting for the pulsar P = 2 and a flux at 1 keV of A = 8.86 x 10 ~^ ph 
cm-2 keV"^ (based on the flux observed by INTEGRAL in 2003; ICotz et allbood ). the 
expected pulsar photon flux is 8.9 x lO^^ x T x E'^ = 5.05 x lO-^E'"^. For LMC X-1 we 
adopt F = 2.5 and A = 0.2 (a typical value in Table 1), for which the corresponding photon 
flux is 0.2 X E^^'^. Equating the two fluxes, the flux of the pulsar dominates over that of 
LMC X-1 only at energies > 1 MeV; even for F(LMC X-l)= 3, the pulsar dominates only 
for energies > 40 keV. 

We now describe a second test performed on our 18 gold/silver spectra which further 
demonstrates that the flux from the nearby pulsar is unimportant. In addition to fitting the 
spectra over our adopted energy range of 2.5-20.0 keV (§§3,4), we also fitted them over the 
ranges 2.5-15.0 keV and 2.5-25.0 keV using precisely the same model and procedures. The 
results are displayed in Figure [9] where it can be seen that the fit parameters and the values 
of xl cire scarcely affected by the choice of the energy range. This is not at all what one 
would expect if the spectrum of LMC X-1 was contaminated by a significant component of 
flux from the pulsar. Rather, one would expect the fit results to vary as the upper energy 
bound is increased from 15 to 25 keV because the pulsar's spectrum is significantly harder 
than that of LMC X-1 (see above). Specifically, one would expect the influence of the pulsar 
to be most apparent at 25 keV. However, the effect of the pulsar on our results, if any, are 
seen to be negligible. In particular, the values of a* (top panel) are completely insensitive 
to the choice of the fitting range. 

As a final direct test, we created a "faked" spectrum of the pulsar in XSPEC using the 
observed flux and power-law index given above and including the off-axis transmission factor 
(T = 0.57). We then subtracted this simulated spectrum from each of the 18 gold/silver 
spectra, redid our standard analysis and found a* = 0.947 ±0.020, which differs by only 0.3 
standard deviations from our adopted result (a* = 0.938 ± 0.020; §5.ip . which was obtained 
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by ignoring the flux of the pulsar. 

We thus conclude, based on the three arguments given above, that the contaminating 
effect of the pulsar on our estimate of the spin parameter of LMC X-1 is negligible (< 1%). 
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Table 1. Preliminary Results for LMC X-1 RXTE Data using phabs(diskbb+powerlaw)^ 



Obs.'' UT MJD'= Tobs Count Rate rms Tin K« Power-Law DPS State*" xl/^oi 

No. (yyyy-mm-(l(i) (s) (cts/s) (10-2)d (kcV) Norm {%) 

(1) (2) (1) (5) (()) (7) (8) (<)) ^l()) (12) (l:5j 
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0, 


.998 


±0, 


.010 


25, 


.28 


± 


2, 


.13 


3, 


,09 


±0. 


,07 


0, 


.34 


±0. 


,06 


53. 


,18 


SPL:TD 





.99/43 


26 


1997-02-08 


50487, 


.468 


9664 


39, 


.49 




0. 


15 


7.58 ±0.4 


0, 


.894 


±0, 


.006 


57, 


.30 


± 


2, 


,92 


3, 


,01 


±0, 


10 


0, 


,23 


±0, 


,05 


66, 


,06 


SPL:TD 


1 


.40/43 


27 


1997-06-18 


50617, 


.157 


5120 


34, 


.28 


± 


0. 


17 


7.86 ± 0.65 


0, 


.873 


±0 


.007 


59, 


.80 


± 


4, 


.50 


3, 


,14 


±0, 


14 


0, 


,23 


±0, 


,07 


68, 


14 


SPL:TD 





.74/43 


28 


1997-07-04 


50633, 


.313 


2464 


31, 


.87 


± 


0. 


,21 


8.42 ± 0.96 


0, 


,846 


±0 


.009 


72, 


.21 


± 


4, 


.21 


2, 


,73 


±0, 


,22 


0, 


.09 


±0, 


,05 


75, 


,25 


SPL:TD 





.76/43 


29 


1997-08-01 


50661, 


.887 


8992 


33, 


.14 


± 


0. 


,14 


7.71 ± 0.42 


0, 


.948 


±0, 


.009 


32, 


.62 


± 


2, 


.75 


3, 


.14 


±0, 


,12 


0, 


.27 


±0, 


.07 


60. 


,74 


SPL;TD 





.99/42 


30 


1997-09-09 


50700, 


.812 


9968 


34. 


.62 


± 


0. 


,14 


7.92 ± 0.39 


0, 


.960 


±0, 


.008 


31, 


.18 


± 


2, 


.41 


3. 


.08 


±0. 


.10 


0, 


.26 


±0. 


.06 


59. 


.73 


SPL:TD 


1 


.24/42 


31 


1997-10-10 


50731, 


.708 


11264 


33. 


.25 


± 


0. 


,13 


9 ± 0.35 


1, 


.012 


±0, 


.010 


21, 


.65 


± 


1, 


.96 


3. 


.12 


±0. 


.10 


0, 


.28 


±0. 


.06 


57. 


.15 


SPL:TD 





.90/42 


32 


1997-11-01 


50753, 


.728 


9008 


33, 


.66 


± 


0. 


,14 


7.63 ± 0.41 


0, 


.989 


±0, 


.008 


29, 


.25 


± 


1, 


.72 


2, 


,86 


±0, 


,13 


0, 


.14 


±0, 


,04 


68. 


,74 


SPL:TD 





.98/42 


33 


1997-11-23 


50775, 


.661 


8336 


35, 


.17 




0. 


14 


8.43 ±0.4 


0, 


.975 


±0 


.009 


30, 


.07 


± 


2, 


.38 


3, 


,05 


±0, 


12 


0, 


,24 


±0, 


,07 


61, 


,82 


SPLiTD 





.99/42 


34 


1998-01-04 


50817, 


.465 


9344 


40, 


.35 


± 


0. 


15 


7.93 ± 0.33 


0, 


.971 


±0, 


.007 


38, 


.55 


± 


2, 


.61 


3, 


,12 


±0, 


,12 


0, 


.27 


±0, 


,07 


66, 


,60 


SPL;TD 


1 


.27/42 


35 


1998-01-25 


50838, 


.750 


9856 


36, 


.03 


± 


0. 


,14 


7.96 ± 0.32 


0, 


.988 


±0, 


.009 


30, 


.17 


± 


2, 


.48 


3. 


.17 


±0. 


.12 


0, 


.28 


±0. 


.08 


63. 


.57 


SPL:TD 


1 


.01/42 



Table 1 — Continued 



vJDS. 


U 1 


i\/r lT\c 


T' . 
obs 


Count Rate 


rms 










i 




Power-Law 


Ur 


oT^ate 


X 


2 /Ar,f 


INO. 


(yyyy-™™-dd) 




(>=) 


(cts/s) 




(iU ) 




yKe V ) 












Norm 




(/o) 












{•^) 




(5) 








(1) 






(J) 






(10) 












36 


1998-02-20 


50864.466 


9728 


37.46 ±0. 


14 


9.44 ± 0.35 


0, 


.994 ± 0.009 


26.72 ± 2.34 


3, 


.12 ± 0, 


,09 


0, 


,33 ±0, 


,07 


56.17 


SPL:TD 


1 


.23/42 


37 


1998-04-07 


50910.312 


9584 


37.91 ±0. 


,15 


8.06 ± 0.36 


0, 


.952 ± 0.009 


34.39 ± 3.21 


3, 


.28 ± 0, 


,09 


0, 


.42 ± 0, 


,09 


55.77 


SPL:TD 





.86/42 


38 


1998-05-28 


50961.227 


9872 


38.37 ±0. 


15 


7.91 ± 0.36 


0, 


.923 ± 0.009 


42.00 ± 3.69 


3, 


.31 ± 0, 


10 


0, 


,44 ± 0, 


10 


56.82 


SPL;TD 





.67/42 


39 


1998-06-28 


50992.190 


9760 


38.25 ±0. 


,15 


7.99 ± 0.34 


0, 


.966 ±0.012 


29.39 ±3.42 


3, 


.42 ± 0, 


,08 


0, 


.59 ±0, 


11 


50.00 


SPL:TD 





.89/42 


40 


1998-08-13 


51038.848 


9200 


38.38 ±0. 


,15 


8.92 ± 0.35 


0, 


.971 ± 0.010 


29.09 ± 2.89 


3, 


.21 ± 0, 


,09 


0, 


.43 ± 0, 


,09 


51.89 


SPL:TD 





.82/42 


41 


1998-09-02 


51058.798 


9264 


35.71 ± 0. 


15 


7.87 ± 0.4 


0, 


.940 ± 0.011 


34.45 ± 3.70 


3, 


.35 ± 0, 


10 


0, 


.45 ± 0. 


,10 


55.07 


SPLiTD 


1 


.01/42 


42 


1998-09-29 


51085.826 


7120 


33.14 ±0. 


,14 


9.33 ± 0.59 


0, 


.964 ±0.013 


25.16 ±3.21 


3, 


.25 ± 0, 


10 


0, 


.41 ±0, 


,09 


49.87 


SPL 


1 


.16/42 


43 


1998-09-30 


51086.525 


2992 


33.84 ±0. 


,20 


10.89 ±0.57 


1, 


.089 ± 0.030 


11.02 ±2.44 


3, 


.12 ± 0, 


10 


0, 


.40 ± 0, 


,09 


40.62 


SPL 





.98/42 


44 


2004-01-08 


53012.266 


8688 


31.24 ±0. 


13 


8.66 ± 0.58 


0, 


.960 ± 0.009 


30.75 ±2.40 


3, 


.06 ± 0, 


12 


0, 


,22 ±0, 


,06 


62.82 


SPL;TD 


1 


.00/37 


45 


2004-01-08 


53012.749 


5360 


30.34 ±0. 


15 


( .OO ± U.OD 


U. 


.y ( o it U.Uli 


zo.DZ it O.U4 


Q 
O 


-t 1 r. 
. ID ± U. 


ID 


0, 


,23 ±0, 


,09 


DO.Ul 


GT3T -T'T^ 


U 


ftft /QV 
.DD/o ( 


46 


2004-01-09 


53013.166 


9264 


31.20 ±0. 


,13 


8.9 ±0.54 


0, 


.977 ±0.010 


25.97 ±2.45 


3, 


.17 ±0, 


,11 


0, 


.29 ±0, 


,07 


58.06 


SPL:TD 


1 


.05/37 


47 


2004-01-09 


53013.801 


3168 


30.84 ±0. 


18 


8.19 ± 0.95 


0, 


.901 ±0.009 


51.84 ±4.36 


3, 


.05 ± 0. 


,26 


0, 


,14 ±0, 


,09 


75.47 


SPL:TD 





.72/37 


48 


2004-01-10 


53014.184 


10160 


30.97 ±0. 


13 


8.61 ±0.65 


0, 


.948 ± 0.008 


33.78 ±2.51 


3, 


.12 ± 0. 


,12 


0, 


.23 ±0. 


,07 


64.73 


SPLiTD 


1 


.15/37 


49 


2004-01-10 


53014.783 


6480 


32.81 ±0. 


,15 


9.2 ±0.59 


0, 


.975 ±0.011 


27.61 ±2.88 


3, 


.17 ±0, 


,12 


0, 


.30 ±0, 


,08 


58.16 


SPLiTD 





.59/37 


50 


2004-01-11 


53015.428 


10800 


31.06 ±0. 


,13 


7.75 ± 0.56 


0, 


.946 ±0.010 


32.07 ±3.10 


3, 


.33 ± 0, 


11 


0, 


.36 ±0, 


,09 


59.00 


SPL:TD 





.82/37 


51 


2004-01-11 


53015.871 


4064 


33.19 ±0. 


17 


11 ±0.61 


1, 


.046 ±0.013 


18.82 ±2.00 


2, 


.83 ± 0. 


,13 


0, 


.18 ±0, 


,06 


57.70 


SPL:TD 





.74/37 


52 


2004-01-12 


53016.340 


10704 


32.75 ±0. 


,13 


10.23 ±0.44 


1, 


.051 ±0.010 


17.18 ±1.45 


2, 


.91 ±0. 


,09 


0, 


.22 ±0, 


,05 


54.11 


SPLiTD 





.85/37 


53 


2004-01-12 


53016.822 


4720 


31.45 ±0. 


,17 


11.74 ±0.59 


1, 


.170 ±0.029 


8.28 ± 1.47 


3, 


.22 ± 0, 


,09 


0, 


.42 ± 0, 


,08 


44.71 


SPL 





.70/37 


^Column density fixed at A^h = 4.6 X 10^^ 


cm ^ . 
































''The first 18 spectra 


, which meet the TD criterion rms 


< 0.075, comprise our favored sample of gold/silver sped: 


ra. 















^Midpoint time of observation. MJD = JD - 2,400,000.5. 

"^Total rms power integrated over 0.1-10 Hz in the power density spectrum (2—30 keV). 

^Normalization constant of the thermal component: K oc Rf^, where Hin is the inner-disk radius. 

^Photon index of the power-law component. 

SDisk fraction defined over the energy range 2—20 keV. 

''See §1 and text in lRemillard &: McClintockl ||2006| ) for state definitions. 
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Table 2: Primary Results for LMC X-1 RXTE Data using PHABS(siMPL(g)KERRBB2) 



No. 


MJD 


a* 


M'° 


f 


r 


fsc 


NDF(%) 




1 -^eatx 


1 


50243.888 


0.933 ± 0.007 


1.256 ± 0.040 


1.557 


2.5 


0.049 ± 0.002 


85.87 


0.89/44 


0.145 


2 


50296.917 


0.884 ± 0.008 


1.679 ± 0.047 


1.551 


2.5 


0.044 ± 0.001 


86.84 


0.56/44 


0.164 


3 


50423.854 


0.962 ± 0.003 


1.183 ±0.021 


1.567 


2.5 


0.035 ± 0.001 


90.29 


1.37/44 


0.158 


4 


50516.454 


0.949 ± 0.005 


1.262 ± 0.032 


1.573 


2.5 


0.038 ± 0.001 


89.36 


1.04/44 


0.157 


5 


50528.108 


0.944 ± 0.005 


1.333 ± 0.031 


1.574 


2.5 


0.032 ± 0.001 


90.93 


0.80/44 


0.162 


6 


50554.291 


0.928 ± 0.005 


1.431 ± 0.032 


1.566 


2.5 


0.037 ± 0.001 


89.35 


0.98/44 


0.161 


7 


50575.107 


0.938 ± 0.004 


1.379 ± 0.031 


1.573 


2.5 


0.036 ± 0.001 


89.84 


0.74/44 


0.162 


8 


50596.944 


0.952 ± 0.006 


1.193 ±0.040 


1.569 


2.5 


0.043 ± 0.002 


87.98 


0.72/44 


0.151 


9 


50597.047 


0.939 ± 0.007 


1.260 ± 0.045 


1.563 


2.5 


0.033 ± 0.002 


90.26 


0.78/44 


0.150 


10 


50638.653 


0.916 ± 0.009 


1.437 ±0.054 


1.550 


2.5 


0.033 ± 0.002 


89.98 


0.66/43 


0.155 


11 


50680.829 


0.945 ± 0.005 


1.401 ± 0.037 


1.572 


2.5 


0.033 ± 0.001 


90.83 


0.75/43 


0.170 


12^* 


50703.106 


0.909 ± 0.007 


1.594 ± 0.043 


1.551 


2.5 


0.032 ± 0.001 


90.44 


0.58/43 


0.168 


13 


50710.736 


0.966 ± 0.005 


1.160 ±0.035 


1.554 


2.5 


0.042 ± 0.001 


88.54 


0.67/43 


0.159 


14 


50794.380 


0.926 ± 0.006 


1.499 ± 0.039 


1.561 


2.5 


0.028 ± 0.001 


91.79 


0.57/43 


0.168 


15 


50884.606 


0.965 ± 0.004 


1.257 ± 0.036 


1.569 


2.5 


0.036 ± 0.001 


90.37 


0.83/43 


0.171 


16 


50939.252 


0.947 ± 0.005 


1.295 ± 0.036 


1.562 


2.5 


0.033 ± 0.001 


90.56 


0.84/43 


0.160 


Yjd 


51014.014 


0.936 ± 0.005 


1.374 ± 0.037 


1.561 


2.5 


0.029 ± 0.001 


91.50 


0.49/43 


0.161 


18^* 


53011.559 


0.940 ± 0.005 


1.273 ± 0.036 


1.551 


2.5 


0.030 ± 0.001 


91.03 


0.76/38 


0.151 


19 




0.938 ± 0.020« 


1.348 ±0.137 


1.562 ±0.008 


2.5 


0.036 ± 0.005 


89.76 ± 1.53 


0.78 ± 0.20 


0.160 ± 0.007 



"Column density fixed at Nu = 4.6 x 10^^ cm"^ and with viscosity parameter a = 0.01. 
''The mass accretion rate is in units of 10^* g s^^. 

^The number of degrees of freedom (dof) in Table 1 is one less because here F is fixed at 2.5. 

''The three strictly TD-state spectra; their average spin parameter is a* — 0.928 ± 0.014. 

^The value of the spin parameter was computed using the one-sided Green's function simpl-1 (§4.2). The 
two-sided Green's function gives an essentially identical result with a mean spin parameter a* = 0.937±0.020. 
If one includes a Gaussian emission hne in the fit to all 18 spectra, phabs(simpL(X'(kerrbb2+GAUSSIAN)), 
the average spin parameter decreases slightly by 0.6 standard deviations to a* = 0.922 ± 0.021 (§5.2). 



Table 3: Comparison of Fit Results for Different Models and Data Samples 



Model 


at 


M 


r 


he 


NDF(%) 


Comments 


I'IIAIiS( SLMI'L'SKEH HUU2)" 


O.I):-!8 ± 0.020 


i.:-!i8 = o.i:-!7 


2.'> 


0.0:«) = 0.005 


89. 7() = 1.5:-! 


18 b])eelia wilii 1 lixcd 


PIIABS(siMPL®KERaBB2) 


0.929 ±0.020 


1.388 ±0.136 


2.76 ± 0.29 


0.050 ±0.016 


88.55 ±2.16 


18 spectra with F free 


phabs(kerrbb2+powerlaw) 


0.953 ±0.018 


1.193 ±0.133 


2.5 




84.11 ±2.28 


18 spectra with F fixed 


phabs(simpli8ikerrbb2) 


0.961 ± 0.033 


1.137 ±0.293 


2.5 


0.051 ±0.018 


86.82 ± 3.62 


53 spectra with F fixed 



"a* = 0.944 ±0.019 for T fixed at 2.3 and a* = 0.932 ±0.021 for T fixed at 2.7. In the extreme case of T = 3.1, 
a, = 0.918 ±0.024. 
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Table 4: Results for LMC X-1 for Five Other X-ray Missions using PHABs(siMPL(g)KERRBB2) 



No. 


Mission & Detector 


a* 


M 


f 


r 


fsc 


NDF(%) 






AfH(1.0E22) 


1 


XMM-Newton EPIC 


0.957 ±0.002 


0.746 ±0.011 


1.618 


2.5 


0.002 ± 0.002 


99.35 


0.96 


0.123 


0.443 ± 0.004 


2 


BeppoSAX MECS 


0.969 ± 0.006 


1.057 ±0.039 


1.628 


2.5 


0.021 ± 0.003 


94.18 


1.03 


0.148 


0.459 ± 0.031 


3 


Ginga LAC 


0.914 ±0.022 


1.487 ±0.136 


1.606 


2.5 


0.023 ±0.005 


93.20 


0.97 


0.159 


0.453 ± 0.093 


4 


BBXRT/ Astro- 1 


0.902 ± 0.028 


0.934 ± 0.067 


1.605 


2.5 


0.099 ±0.019 


71.90 


1.06 


0.096 


0.575 ±0.015 


5 


ASCA CIS" 


0.778 ± 0.012 


2.287 ± 0.048 


1.643 


2.5 


0.079 ± 0.005 


78.64 


0.70 


0.176 


0.512 ± 0.006 



"In order to fit this Ginga LAC spectrum, it was necessary to add a Gaussian emission-line component, 
PHABS(simpL(X'(kerrbb2 + GAUSSIAN)); the central energy of the Hnc, its full-width at half maximum and 
equivalent width arc respectively 7.0 ±1.5 keV, 1.3 ± 0.6 kcV and 0.51 keV. 

''ASCA SIS gives a similarly low value of spin a* = 0.742 ± 0.012. 
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Fig. 1. — POWERLAW versus SiMPL. The left pair of panels show unfolded spectral fits to 
a PCA observation of LMC X-1 (No. 5 in Table 1) and the pair on the right show the cor- 
responding unabsorbed models. All of the fits were performed using DISKBB in conjunction 
with POWERLAW (top panels) and SiMPL (lower panels) . The composite model is represented 
by a solid line, the emergent disk component by a dotted line, and the Comptonization com- 
ponent by a dashed line. Note the strikingly difi^erent behaviors of this latter component in 
the unabsorbed models. 
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Fig. 2. — Spin parameter versus observing time. The dashed hne indicates the mean value of 
the spin parameter for the 18 gold/silver spectra. As indicated, the scatter about the mean 
- an effective measure of the observational uncertainty - is ~ 2%. The small error bars 
shown reflect only the counting-statistical errors associated with the many counts detected 
per observation (> 10^ counts; Table 1). 
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Fig. 3. — Histograms of the fitting parameters for our adopted model 
phabs(simpl(8)KERRBB2) with r fixed at 2.5. The sohd histograms are for the 18 
gold/silver spectra and the open plus sohd histograms are for all 53 gold/silver/bronze 
spectra. 
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Fig. 4. — Correlations between pairs of fitting parameters for our adopted model 
phabs(simpl(8)KERRBB2) with r fixed at 2.5. Data for the 18 gold/silver observations 
are plotted as filled circles and for the 35 bronze observations as open circles. The mass 
accretion rate M is in units of 10^^ g s~^, and rj is the efficiency for converting rest mass into 
radiant energy. 
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Fig. 5. — The scattered fraction /gc (open circles, scale on the left) and the naked disk 
fraction NDF (filled circles, scale on the right) versus the rms power (0.1-10 Hz) for our 
standard model PHABS(siMPL(g)KERRBB2) with V fixed at 2.5. The solid vertical line sepa- 
rates the data for our adopted sample of 18 gold/silver spectra (rms power < 0.075) from 
that of the 35 bronze spectra. We find very similar correlations for F = 2.3 and F = 2.7, i.e., 
the scattered fraction generally increases with increasing F while the disk fraction decreases. 
This has the effect of shifting the plots up or down relative to the F = 2.5 plot shown here. 
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Fig. 6. — Exploration of the relationship of the classical TD rms criterion (rms power 
< 0.075) to the scattered fraction. The spectra are plotted in order of increasing rms power. 
Top: The first 18 data points with rms power < 0.075 (indicated by the vertical dashed 
line and solid horizontal line) correspond to our gold/silver spectra and the remaining 35 are 
the bronze spectra. Bottom: The horizontal dashed line indicates the scattered-fraction 
cut line at 4% that corresponds to the rms power = 0.075 line (see text). The spectra with 
fsc < 4% are indicated by filled circles in the top panel. As shown, there is a very good 
correspondence between favored spectra with low rms power and those with low scattered 
fraction. 
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Fig. 7. — Effect on the spin parameter a* of varying the input parameters M, i, and D for the 
viscosity parameter a = 0.01. (a) The upper panel shows a normal distribution for the black 
hole mass and the lower panel shows a* versus mass M for 3000 sets of parameters drawn 
at random. The central filled circle indicates our estimate of the spin (d^o = 0.937lo;og°) 
obtained from these simulations. (b,c) Same as panel a except now for the parameters of 
inclination angle and distance, respectively. That the uncertainty in the mass is the dominant 
source of observational error is indicated in panel a by the significant correlation between 
a* and M. (d) Histogram of spin displacements for 3000 parameter sets. The vertical solid 
line indicates the median value of the spin in the simulated sample, a*o = 0.937. The two 
dashed lines enclose 68.27% (Icr) of the spin values centered on the solid line and imply an 
observational uncertainty of (+0.040, -0.061). 
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Fig. 8. — Combined error analysis that considers both of our fiducial values for the viscosity 
parameter. The thin solid fine is for a = 0.01 and the dashed fine is for a = 0.1. The sum of 
these two smaller histograms forms the large histogram (thick solid line), which shows the 
spin displacements for a total of 6000 parameter sets. The vertical sohd fine indicates the 
median value of the spin determined by these simulations: a^,o = 0.92. The two dashed lines 
enclose 68.27% (Icr) of the spin values centered on the solid line and imply an observational 
and modcl-paramctcr uncertainty of (+0.05, -0.07). The two dotted lines enclose 95.45% 
(2(j) of the spin values and imply an uncertainty of (+0.06, -0.18). 
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Fig. 9. — Evidence that the nearby pulsar PSR B0540-69 does not significantly contami- 
nate the 18 gold/silver spectra (see Appendix B). In making this test, we used the model 
PHABS(simpL(S)KERRBB2) and allowed the power-law index F to be a free parameter. The 
top four panels show our standard fitting parameters that are defined in §4.3 and the follow- 
ing panel shows the naked disk fraction NDF, which is defined in §4.2. The bottom panel 
gives the reduced chi-square of the fit. The fits were performed over three energy ranges: 
2.5-15.0 keV (open circle), 2.5-20.0 keV (plus symbol) and 2.5-25.0 keV (asterisk). The 
mean value of each parameter for our nominal energy interval of 2.5-20.0 keV is indicated 
by a horizontal dashed line and is given directly above the dashed line inside each panel; the 
mean values for the other two energy intervals are indicated by the symbols plotted at the 
far right of the figure. The vertical dashed lines indicate the three gold TD spectra (Nos. 
12, 17 & 18). 



